Ion acceleration in "dragging field" of a light-pressure-driven piston 
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We propose a new acceleration scheme that combines shock wave acceleration (SWA) and light 
pressure acceleration (LPA). When a thin foil driven by light pressure of an ultra-intense laser 
pulse propagates in underdense background plasma, it serves as a shock-like piston, trapping and 
reflecting background protons to ultra-high energies. Unlike in SWA, the piston velocity is not 
limited by the Mach number and can be highly relativistic. Background protons can be trapped 
and reflected forward by the enormous "dragging field" potential behind the piston which is not 
employed in LPA. Our one- and two-dimensional particle-in-cell simulations and analytical model 
both show that proton energies of several tens to hundreds of GeV can be obtained, while the 
achievable energy in simple LPA is below 10 GeV. 

PACS numbers: 52.38.r, , 42.65.Re, 52.27.Ny, 52.65.Rr 



Laser-driven ion acceleration attracts more and more 
attention nowadays because of its potential to realize 
table-top size ion accelerators, which would greatly re- 
duce the expense and occupying area compared to con- 
ventional accelerators. In most proposals, intense/ultra- 
intense laser pulses are employed to stimulate strong elec- 
trostatic fields with amplitude several magnitudes higher 
than that of conventional methods. These fields move 
at high velocities thus ions can be continually acceler- 
ated to large energy in a short distance. In Target Nor- 
mal Sheath Acceleration (TNSA) ^^H], the sheath field, 
formed at the target back by laser heated super-thermal 
electrons, moves at several to tens of the ion acoustic 
speed. Ions co-move with the field and are accelerated to 
tens of MeV, while energy above hundred MeV is quite 
a challenge because of the weak scaling law versus laser 
intensity. In Light Pressure Acceleration (LPA) [9Hr6j, 
electrons in the thin plasma foil are pushed inward by 
the laser pressure and build up an intense charge separa- 
tion field, by which protons are accelerated. As the foil 
is driven as a whole, protons move along with the elec- 
trostatic field and are successively accelerated to energies 
that theoretically may reach GeV. Nevertheless, even in 
simulations it is difficult to overcome the ten GeV bar- 
rier for protons with LPA, since the accelerating gradi- 
ent drops quickly due to the relativistic Doppler effect. 
As soon as the protons become relativistic, their energy 
grows very slowly with time (~ t^/^). 

Recently, a breakthrough in collisionless shock wave ac- 
celeration [17- 22J has been reported [19j. It suggests that 
a multi-pulsed picosecond (ps) laser pulse can efficiently 
excite a collisionless shock and maintain its velocity in a 
density-decaying plasma target. The produced protons 
have much larger energy and lower energy spread than 
predicted by the hole-boring shock [20j . When the shock 
propagates like a piston, some protons are trapped and 
refiected to a speed determined by the energy density gra- 
dient. The exponentially decaying density profile of the 



target is important for keeping the shock velocity high 
and constant. However, plasma itself also expands at the 
ion acoustic speed simultaneously. It seems inevitable 
that the pressure gradient will be further weakened at a 
later stage. This would put a limit on the final energy 
gain. From other side, the shock wave velocity is limited 
by the Mach number, which cannot be very high [21j. 

Generally, the shock front may be considered as a fast- 
moving piston carrying a strong electrostatic field. Imag- 
ine, one could free the limitation on its velocity, e.g., the 
shock becomes relativistic [23j. Then, the trapped and 
refiected protons would reach a relativistic 7— factor of 
about 7 = 27^, where 7^ is the relativistic factor of the 
shock wave. This is a way beyond the simple thermal 
shock wave acceleration. An intense electrostatic field 
moving at relativistic velocity is required. One possibil- 
ity might be the plasma wakefield. It is works very good 
for electron acceleration, however it may be not intense 
enough to trap protons, although there are some propos- 
als to trap pre-accelerated protons and accelerate them 
further by the wakefield to very high energes [24j. 

We propose here that the electrostatic field of LPA can 
serve as the perfect relativistic piston. Its velocity is close 
to the light speed and the charge separation field is suf- 
ficiently high to trap background protons. In this paper, 
we develope an analytical model and use particle-in-cell 
(PIC) simulations to show that when a thin foil driven by 
ultra-intense circularly polarized (CP) laser pulses, with 
peak amplitude ao = eEL/rrieUJoc from 50 to 200, passes 
through an underdense plasma region, refiected protons 
can achieve energes up to hundreds of GeV. Here e and 
rrie are fundamental charge and electron mass. El and 
uq are electric field amplitude and angular frequency of 
the laser pulse, c is the light speed, respectively. 

The expression of the accelerating field amplitude as 
a function of time is derived, which, together with the 
foil momentum equation, well describes the movement of 
the background protons and the trapping condition. The 
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Figure 1: The electrostatic field driven by a relativistic CP 
laser propagates in the underdense background plasma, traps 
background protons and reflects them to ultra-high energy. 



scaling laws indicate that the peak energy increases with 
t^/^ and ttg. A proof-of-principle two-dimensional (2D) 
simulation is also performed to show the validity of our 
analysis in a multi-dimensional geometry. 

The mechanism is sketched in Fig. [T] . As known in 
LPA, the laser-driven foil is accompanied by an intense 
charge separation field. When it propagates through the 
low density plasma, background protons can be contin- 
uously trapped in certain conditions. The foil plays the 
role of a relativistic piston. In the frame of the relativis- 
tic piston, background protons clash towards the electro- 
static field at the foil velocity, then are slowed down by 
the field in and behind the foil and finally reflected. 

In LPA, protons are accelerated by the field local- 
ized in the ultra-thin skin layer only, leaving the enor- 
mous "dragging fiel'" (DF) unused. In our scheme, the 
background protons pass through the skin layer and are 
greatly accelerated by the DF, which contains much more 
electrostatic potential than the field inside the layer. As 
a result, the reflected protons can obtain energy a magni- 
tude higher than foil protons in the simple LPA regime. 

To demonstrate the mechanism a ID PIC simulation 
is performed by the VLPL code [25j. A dense thin foil, 
with density of n/ = IbOric and thickness of df = 0.57Ao, 
is driven by a CP laser pulse with peak amplitude ao = 
150, where Uc = rrieOjQ / Aire'^ is the critical density and 
Ao = 0.8/im is the laser wavelength. For simplicity, we 
used a trapezoidal laser pulse . Its amplitude increased 
linearly to the maximum in 2To and then stayed con- 
stant (To is the laser period). Here the foil thickness is 
optimized for LPA [26j . A background plasma with den- 
sity of rib = O.OOlnc is located behind the foil. Fig. |2|a) 
shows distributions of the electrostatic field and momen- 
tum of background protons. The field decays quickly at 
first and much more slowly later. It is so intense, well 
beyond Eq = rrieUJoc/e even at the end of the interac- 
tion, that background protons are gradually trapped and 
reflected to about 70 GeV in 3mm, while protons in the 
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Figure 2: a) Distributions of electrostatic field (blue solid) 
and momentum of background (BG) protons (red dotted) at 
different stages with ao = 150, n/ = 150nc, dj — 0.57Ao and 
Uh — O.OOlric. b) Distributions of laser field (black solid), 
electron density (red solid) and proton density (blue solid) 
at t = 500To. c) Evolution of the peak electrostatic fields of 
ao = 50, 100 and 150, from analytical model (blue solid) and 
one-dimensional simulations, d) Evolution of the maximum 
energy of background protons with ao = 50, 100 and 150, 
from analytical model (blue solid) and simulations. Here the 
dashed lines show that the peak energy increases with ^ t 
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front foil have the peak energy below 10 GeV. 

Trapping and accelerating of charged particles are de- 
termined by the velocity, amplitude and length scale of 
the DF structure. Since the DF co-moves with the foil, 
its dynamics can be derived by the momentum equation 
of the piston [llj 



dt mi NfDf 1 + /3/' 



(1) 



here /3/ is the foil velocity normalized by c, 7/ = (1 — 
^2j-i/2^ and rui is the ion mass. Nf and Df are initial 
density and thickness of the foil normalized by the cirit- 
ical density ric and laser wavelength Ao- All lengths and 
the time are normalized by Ao and To, respectively. The 
velocity of the DF structure is thus Pe =Pf^ ^ind can be 
obtained by solving Eq. ([T]). 

More critical is the evolution of field amplitude. For- 
mer researches suggested that considering the balance of 
the electron thin layer between the light pressure and 
electrostatic force, the peak amplitude of the DF should 
be proportional to the Doppler factor (1 — /3/)/(l + /3/) 
[27J. This means that it would drop by l/47j in the 
relativistic regime. In the above simulation 7/ !^ 5 at 
t = 5OOT0, and the DF amplitude should decrease by 
some factor 100. However, as seen in Fig. [2|a) The DF 
at t = 5OOT0 is still above E/Eq > 15, roughly one tenth 
of initial maximum field. A new scaling should then be 
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deduced. In Fig. [2jb) the distributions of laser field and 
electrons at t = SOOTq are presented. The laser pondero- 
motive force acts on only a small portion of the electrons 
at the surface and not on the whole electron layer. This 
is because the light pressure is severely weakened when 
the foil becomes relativistic so that it cannot confine all 
foil electrons. The situation is then similar to the elec- 
trostatic shock or the hole-boring process [9j, except the 
whole system is moving fast. 

Since the foil velocity varies slowly with time, it is rea- 
sonable to assume the piston frame as an inertial ref- 
erence system, where relationships in Ref. [9j can be 
applied after Lorentz transforming all quantities. Con- 
sidering the balance of the electron skin layer instead of 

the whole electron layer, the peak DF is E^q ~ \[^lad 
, where P^^^ is the light pressure in the moving frame 
K . The longitudinal electric field and light pressure 
turn out to be Lorentz invariants, thus the peak am- 
plitude in laboratory frame is conveniently obtained by 
£^xO - VP^d - + On the other hand, 

the ponderomotive force acting on electrons is approxi- 
mately |e/3g X B^l ^ eE]^ , therefore E^q ~ E^^ giving 
the normalized peak amplitude of 




(2) 



The both estimations above show that the field scales 
proportional to the square root of the Doppler factor, i.e., 
about 1/27/ . This is a much weaker decay than formerly 
assumed. We compare the results from Eqs. ([T]) and ([2|h 
with ID PIC simulation results in Fig. |2j One can see 
that for a=50, 100 and 150, Eq. [2] describes the simula- 
tion results properly. This new scaling is very important 
to guarantee the DF structure can be maintained for a 
long distance and does not vanish in relativistic regime. 

The dynamics of background plasma protons encoun- 
tered in the field of the piston are then derived by 



dt 



rur, 



(3) 



As seen from Fig. [T]and Fig. [2|a), the DF drops almost 
linearly with the distance shortly behind the foil. The long 
weak tail far behind the foil is neglected since protons 
located there are no longer trapped. We simplify the de- 
scription, by assuming the charge separation field behind 
the skin layer decays to zero in a length scale of 
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elsewhere 



(4) 

Here the X5, /35, 75 are position, normalized velocity and 
7— factor of the background protons. The scale length of 




Figure 3: Scaling law of the maximum achievable momen- 
tum versus peak laser amplitude from analytical model (black 
square) and simulations (blue pentacle). The fitted red line 
indicates the power-law of ^ ao- 



the electrostatic field localized inside the skin layer is so 
small that we ignore it. Simulations suggest that the DF 
scale is about dE ^ 40/im. Protons initially located at 
different positions of x^lt^o, which are also the injected 
positions, obtain final energy when x^ = Xf. So the 
maximum energy of all reflected protons as a function of 
time can be analytically derived by changing the X})\t=o 
in Eqs. d3| and Q. 

Eqs. (jlf-Q offer a complete description of the mecha- 
nism. Peak momenta of the trapped protons at different 
simulation times and from analytical model are shown in 
Fig. [2|d), for a=50, 100 and 150, respectively. Again, 
the analytical results and the simulations are in a good 
agreement. Simulation data is also fitted, showing that 
the peak momentum increases with time according to the 
power law cx t^/^ , much faster than in the simple LPA 
(oc t^/^) [llj. The great difference is due to the fact 
that the light pressure - that drives the foil - decays as 
1 / 47 J while the DF accelerating the background protons 
decays much slower, as 1/2^ f . Within a few millime- 
ters the protons are accelerated close to lOOGeV, while 
in LPA it is nearly impossible to produce protons beyond 
ten GeV even in ID simulations. 

In DF acceleration, the piston is continuously accel- 
erated. Thus, the later local protons are injected, the 
higher energy they obtain. However, the DF amplitude 
also decays meanwhile so that after certain point it be- 
comes not intense enough to trap stationary protons. The 
trapping condition is obtained by balancing kinetic en- 
ergy of the incoming protons and the DF electrostatic 
potential in the relativistic piston frame 



-eE^Q{t)dE = (7/ - 



(5) 



where E^Q{t) = Exo{t) and d^ = jfdE is the peak ampli- 
tude and length scale of the DG field in the moving frame. 
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The critical foil 7— factor beyond which background pro- 
tons are no longer trapped is then 



7/c 



1 



irrrie a^DE 



(6) 



For ao = 150 and De = dE/^o = 40, the critical value is 
about 7/c ^ 6 . After the foilaAZs 7— factor reaches the 
critical value, latish protons would pass through the DF 
structure without being trapped any more. This defines 
the maximum energy obtainable by this mechanism. In 
Fig. |3]the scaling law between the maximum obtainable 
energy and ao from simulations is compared with the 
model. The simulation results are only shown for ao = 50 
and 100. For higher amplitudes, it would take too much 
time to reach the final stage (about 10^ lijm). Fig. [s] 
indicates that the maximum energy is proportional to 
the square of the laser amplitude ~ aQ. This mechanism 
allows to accelerate protons to hundreds of GeV. 

In ID simulations we chose very low background den- 
sities to ensure that the DF is not affected. In the multi- 
dimensional geometry, the situation becomes more com- 
plicated. Firstly, instabilities may destroy the thin foil 
during the interaction. This could be restrained by using 
a highly relativistic laser pulse, e.g., ao > 100 in our case. 
Secondly, the foil deformation due to the laser intensity 
distribution will cause diffraction. Refiected parts of the 
laser fields interfere with each other behind the moving 
piston. The DF structure is well formed at the begin- 
ning, however when the perturbation arrives at the axis, 
it may breakdown before any protons being trapped. 

To address this, we increase the background plasma 
density. Background electrons are snow-plowed by the 
driving laser and pile up inside the piston. As a re- 
sult, higher background density could induce two effects. 
Firstly, the relativistic piston is slowed down. Secondly, 
the DF might be enhanced by the additional charge sep- 
arating field between piled-up electrons and left-behind 
protons. The higher background density easies the trap- 
ping of local protons. The 2D simulation results are shown 
in Fig. [4j The laser pulse has the same trapezoidal pro- 
file as in ID simulations in time domain and a super- 
Gaussian distribution transversely ~e~^^/^^^ , where 
Wy = 35/im. To restrain instabilities, the laser peak am- 
plitude is increased to ao = 200. The diffraction and 
refiection of the incident pulse is minimized by using 
density matching foil, i.e., the foil is with density dis- 
tribution of nf{y) = riQe'^y^'^y^ [16j. The foil density 
and thickness are Uf = 80nc and df = 0.62Ao while the 
background plasma density is increased to 77,5 = O.lnc- 

Fig. |4ja) shows distributions of the DF, electrons in 
the foil and background protons at t = 6OT0. It is seen 
that the DF is well formed behind the thin layer. Before 
the instabilities could develop perturbing the accelerating 
field, some of the background protons are trapped and ac- 
celerated. As seen in Fig. |4]^b), the maximum energy of 
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Figure 4: Two-dimensional results with ao = 200, Uf — SOric, 
df — 0.62Ao and Uh — Clric. a) Space distributions of the 
foil protons (black dotted), background protons (red dotted) 
and the electrostatic field at t = 6OT0; b) Momentum dis- 
tributions of the foil protons (black dotted) and background 
protons (red dotted). The right panel is the energy spectrum 
of background protons. 



refiected protons at t = 500To is about 25 GeV, approx- 
imately 4 times the peak proton energy in LPA. Some 
foil protons left behind the layer are also trapped and 
accelerated to high energy. The total energy of refiected 
background protons is about 3% of all foil protons, which 
is quite considerable concerning the low background den- 
sity. One would expect higher energy efficiency by using 
denser background plasma. 

In conclusion, a method to laser accelerate protons be- 
yond tens of GeV is proposed. It takes advantage of the 
electrostatic field dragging behind the fiying foil in LPA 
to trap background protons. The dragging field keeps 
trapping and accelerating protons to energies greatly 
higher than that obtained in LPA. The analytical model 
shows that the electrostatic field decays more slowly than 
formerly assumed, resulting in larger energy increasing 
rate (- t^^^).T\ie final maximum energy scales as ~ aQ. 
2D simulations proved the proposal, generating protons 
with peak energy of about 25GeV. 
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